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ABSTi^CT 


The present day EHV/UKV pov/er system network requires 
highly reliable/ sensitive/ selective and fast protective 
schemes. The overcurrent relays were replaced gradually by 
distance relays for transmission line protection due to 
their inherent drawbacks. The static relays using solid 
state cortponents were installed in sixties. The first idea 
of enploying digital computer for the protection of pov/er 
system networks was proposed during 1969. Till today/ a 
number of corrputer oriented real-time protection algorithms 
have been proposed. In this thesiS/ the overview of the 
digital computor distance relaying algorithms developed 
uptill now has been presented in brief. The predictive 
calculation of the impedance from the relay location to the 
point of fault is presented and tested for several faults/ 
a good accuracy has been obtained in the impedance calculation. 
Different types of protective schemes have been realized 
and tested by digital simulation. 
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CHAPTER 1 


INTRODU CTION 


Power system network is- becoming more and more 
conplex today thus requiring the protection schemes that are 
highly reliable, sensitive, selective and also, fast in 
operation. A transmission line which is the link betv/een 
the source and the load at E?IV/UHV level requires special 
attention in the field of protection. In the past, over- 
current relays were generally used for the protection of the 
transmission lines. But, these relays are being replaced 
gradually by distance relays due to certain demerits such as, 
change of balance point with the type of faults, cnange in 
generation capacity and, switching operations etc. Of the 
several electromechanical measuring units available, the 
induction cup unit proved to be the best, because of several 
advantages. However, with the advent of semiconductors, the 
static relays using solid state corrponents were installed in 
sixties. The static relays possess several meirtis such as, 
greater sensitivity, lower burden, fast . in operation, absence 
of contacts and mechanical motion and hence less maintenance 
and immunity from 'v^ib-ratA'^® and shocks due to external causes 
such as movement, electric locomotives, earthquakes etc. The 
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static relays are being used increasingly in recent yearS/. 
particularly for the protection of EHV/UHV systems where 
increased sensitivity^ reliability and speed are iuportant. 

Use of ICs in place of transistors and discrete corrponents 
leads to corrpact, simpler, more reliable and economical 
relay ■unit. Studies have sho'wn that a given power system 
can transmit larger amount of power without loss of synchro- 
nism, i.e, stability when a fault is cleared in a shorter time. 
With the advent of fast operating rel and fast acting CBs, 
now, it is possible to remove a fa'ult on a transmission line 
within two cycles (i.e. 40 ms). Higher speed clearing of 
faults also minimizes the possibility of permanent damage to 
the conponents, fire and h\iraan risks etc. 

Digital computers have been applied to both off line 
and on-line ■ttaato pertaining to power system studies (load flow, 
short circuit and stability), control and operation. However 
■their use for the protection of the power system networks is 
of recent origin, the first idea was proposed during 1969 in 
this field. Till today, a member of computer oriented real-time 
protection algorithms have been proposed. Extraction of the 
fundamental conpon^nts from the complex post— fa\ilt waveforms 
with the aid of analog and/or digital filters is necessary wi^th 
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■the abo'VQ algori-thms. The digital computer protection of 
the transmission lines possesses several mertis such as 
flexihili'ty^ i. e, any desired threshold characteristic can be 
obtained ■with the same hardware or wi'th the minimum changes 
in the hardware, self checking properties, adaptive capabili^ty 
and data in -ter face access facilities. 

In this thesis, the overview of the digi^tal conputer 
distance relaying algorithms developed uptill now has been 
presented in brief. ©Ut of the several algorithms proposed, one 
of them is presented in detail and also tested on DEC 1090 at 
IIT Kanpur for different possible faul^fcs and all possible 
threshold characteristics have been developed. 

The use of digital corrputers and microprocessors for 
protective relaying purposes has been drawing the attention of 
researchers and protection engineers since late 1960* s. 

G.D, Rocbefeller Cl 3 in 1969, proposed a conple^te groups of 
algorithms for the protection of "the transmission line and 
other power system conponents using digi-tal coirputers. The 
algorithms proposed up till now can be classified into two 
categories, distance relay algori^thms and travelling-wave relay 
algoritlvns. The distance relay algorithms involve the deter- 
mination of irtpedance, from the relay location -to the point 
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of fault, from the fundamental components of voltage and 
current v/hich are extracted from conplex post-fault waveforms 
by analog and/or digital filters. The travelling-wave relay 
algorithms utilize the complex post-fault waveforms as they 
are, for making relaying decisions by means of travelling wave 
techniques. Slemon et al [2] and Ramamoorty [3 ] described the 
determination of the f-ondamental conponsnts in phasor form from 
an ensemble of samples collected over one full cycle by using 
Fourier Analysis technique. Mann and Morrison [ 4 ] described 
the predictive calculation of impedance and its phase angle 
from a much fewer number of samples of voltage and current and 
their tirae derivatives, in another paper [ 5 ] , they proposed 
a relaying of a three phase line using this technique, 

P, Gokul et al [6] described the predictive calculation of 
impedance and its phase angle and also quadrilateral character- 
istic was realised. The predictive technique was presented 
also by Gil crest et al [ 7 ] who proposed the use of first and 
second time derivatives of the sarrpl es, and, by Gilbert and 
Shovlin [ 8] who proposed the use of samples and sampling interval. 
Cars and Jackson [ 9 ] described the use of two digital orthogonal 
notch filters with sine characteristics to determine the magni- 
tude and phase angle of the fundanvental components from samples 
taken at four equally spaced time intervals over the fmdamental 
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period, and, the use of an analog low pass RC filter with a 
cut off frequency of 85 Hz to band limit the voltage and 
current signals prior to sampling. Hope and Umamaheswaran[lO ] 
described the use of odd and even square waves, besides sine 
cosine functions, for the extraction of fxandamental frequency 
components, Hope et al [11] and Thirupathaiah et al [ 12 ] 

f 

have also enployed Fourier Analysis techniques with one cycle 
data window, whilst Phadbe et al[ 13,14], and vJisznicwski [l5 ] 
and, Sunnak et al [ 16 ] have employed these techniques vrith 
half a cycle data window. Some authors proposed the transmission 
line modelling by a series R-L circuit, the numerical solution 
of the differential equations provides the value of R and L, 
Njttlnnes and Morrison [ 17 ] proposed the integration of the deffe- 
rential equations over tv/o successive sampling intervals to 
generate adequate number of equations to solve for R and L, and 
used the trapezoidal rule to evaluate the integrals. Poncelet 
[ 18 ] / and Bernard and Bastide [ 19 ] treated the deficiency in 
this modelling as an error and solved for R and L subject to 
the minimisation of the sim of squares of these errors over a 
certain niomber of successive sampling intervals. Smolinski 
[ 20 ] modelled the transmission line by a single Pi section- 
including the shunt parameters, and solved the resultant diffe- 
rential equations by replacing them by finite differences and 
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using four saxrpla sets to obtain the value of R and L, In 
the above modelling^ the filtering is not needed as the cause 
of h,f, transients, shunt capacitors, has been included in 
the modelling. Ranjbar and Oary[ 21 ] presented a digital 
harmonic filter in which the integration is carried out over 
a certain manaber of overlapping sxibintervals with preset limit 
while the differential equations are being solved by numerical 
integration method, Sachdev and Baribean [ 22 ] assumed the 
post fault waveform to be comprising of a decay D.C. offset, 
the fundamental and a certain nvimbers haritionicc and determined 
the fundamental component by applying the least square error 
method, Johns and Martin [23] presented a finite transform 
method in which the fundamental frequency impedance can be 
determined by carrying out the filtering process in the frequency 
domain, Girgis and Brown [24, 25] developed two state and 
three state Kalman filters to extract the fundamental voltage 
and current phasor respectively, and reported, that, the error 
involved is less than one percent after half a cycle. Sanderson 
and Wright [ 26] have dealt v;ith series compensated lines. They 
modelled the line bya series RLC circuit and solved the resxxl- 
tant differential equations by integrating them over three 
successive sampling intervals. 
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Vitin's correlation method based on travelling wave 
theory [ 27 ] is suitable for computer application also. Tahagi 
and Yamakosi [ 28 ] described the microcomputer application of 
the differential current relay [ 29 ] v/ho se operation is based 
on the travelling wave theory, Deshilcachar and L.P, Singh[30] 
developed the relaying algorithm based on travelling wave 
phenomena/ which is very fast in operaticn . 

NoW/ a chapterwise description of the work carried 
out and reported in this thesis is given belov7. 

Chapter two presents the distance relay algorithms 
proposed by several authors and also the theory and principles 
of filters to extract the ftmdamoital components from the po st 
fault waveforms which, in general, are not sinusoidals as 
they contain in addition to the fundamental coirponents, d.c, 
offset and transients corrponents also. 

Chapter three deals with the digital simulation of a 
proposed three phase relaying scheme with all possible 
threshold characteristics, 

Finally, chapter foxir concludes vmth the findings 
of the present worlc and scope for future extension. 



CHAPTER 2 


DIGI TAL SOIEMES ^_Og_Ti^a^ISS ION NETWORKS 

2.1 INTRODUCTION 

The need, for fast and reliable protective relaying 
schemes for transmission lines, especially at SHV/UHV voltage 
level was felt a long back. In fact, the installation of 
minicomputers and microprocessors for the purpose of protec- 
tion of power system networks, which is a step in this 
direction, is going to be reality in near future. Digital 
protection of transmission lines possesses several merits such 
as faster relay operation and flexibility to obtain any type 
of desired threshold characteristics with the same hardware 
or with the minimum changes in hardware and also self checking 
against any failure in hardware. 

The work in the field of protection of power system 
networks using digital coirputer was started in late 1960s and 
since then, a nxrriber of algorithms have been proposed for the 
protection of transmission lines. The algorithms proposed so 
far can be classified into two groups; 

1. distance relay algorithm based on fundamental voltage and 
currant, 

2, travelling wave relay algorithm based on travelling wave 
phenomenon. 
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Both these algorithms need appropriate type of filter and 
filtering circuit which is being discussed in the following 
paragraph, 

2 . 2 PILT3RING 

The post fault waveforms, during the first one or 
tvro cycles following the occurance of a fault, corrprise •fa 
power frequency fundamental corrponent along with an ejqso- 
nentially decaying d.c, component as well as high frequency 
transients. In addition, if the line fault is near generating 
station, the post fault line voltage and current waveforms 
will also contain subtransient and transient power frequency 
components. The distance relaying algorithms involve in 
computing the power frequency impedance between relay location 
and the point of fault from the fundamental coirponents of 
voltage and current obtained by suppressing other corrponents 
through analog and/or digital filters and line modelling. 

This filtering introduces a time delay in relay operation. 
Analog and digital filters, which have been proposed and/or 
being used for this purpose, are described in, the following 
sections* 

2,2,1 ^al og F ilters - [4] 

Iron cored current transformers with resistive burdens 
transmit the full amount of offset in the primary f atilt current. 
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This conponent often takes inore than 5 cycles to decay on the 
present systems with high X/R ratios. The voltage can also 
have an offset conponent, but this is considerable smaller 
than that on the current waveform. The relevant expressions 
are well known [32] and are repeated here for convenience, 

-R- t/L 

= ^Ipk SinCwt-i-c-qJ^)- e Sin(a-92_) (2,1) 


= ^lok SinCwt + a _ 


c, + <p 

1 


Sin (a-^^) sin -R^t/L^ 


Sin P, 


( 2 . 2 ) 


In these equations, a is the phase of fault incidence, the 
subscript 1 refers to primary (system) quantities, and L» 
refers to line quantities. 


If an ideal C.T, is connected to a secondary burden 
having the same X/R ratio as the primary circuit, then the 
voltage across the burden will be purely sinusoidal. 


For a C.T, of secondary burden R 2 and L 2 , 


V 


^2 

- ^ V2 ^2 dF*^ 


(2.3) 


where subscript 2 refers to CT secondary 
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Prom above equaticsi, the coefficient of the exponential term 
is zero if. 



Thus, dc offset in the current signal can be eliminated using 
mimic impedance in the CT secondary, to match primary -X/R 
ratio to that of the secondary. But, the exact cancellation 
of the dc offset is not possible for all faults. The primary 
X/'R. ratio to be matched is that of source plus line to the 
fault point, and since in general, the source X/'R is not equal 
to the line X/R, the overall primary X/R is a variable quantity, 
dependent on, how far, along the line the fault occurs. 

RC low pass filters, with a suitable cut off frequency, 
have been used to eliminate the high frequency conponents in 
conjuncticn with some of the algorithms The value of 

the highest frequency required is of considerable inportance 
since, it is necessary for 

f > 2 f 

s max 

where, £_ ® sanpling frequency and 

f =s maxiroTom frequency present in sanpled signal. 
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The cut off or half power frecjuency f * for a sirrple RC low 
pass filter/ is given by. 



The 10—90% rise time tr of this filter to a step input is 
given by, 

tr = RC (- In 0.1 + In 0.9) 

= 2.2 RC (2,6) 

Therefore, by combining above tv/o equations (i.e, eqns. 2,5 
and 2.6), we get, 

f ^ (2.7) 

tr 

The minimum rise time required can be determined from an 
examination of the protection speed requirements. 

2,2,2 Digit al F ilters 

The digital filter is a digital system that can be 
used to filter discrete time signals. It can be inplemented 
by means of software (computer programs) or by means of dedi- 
cated hardware, and, in either case, it can be used to filter 
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real-time signals or non-real time signals. The various 
types of digital filters^ which have been proposed so far, 
are described below. 


2. 2, 2,1 Notch Filters [s] 


The band of frequencies around 50 Hz will be passed 
by a filter with a 50 Hz sine wave iirpuLse response. The 
response must be finite in length for obvious reasons and 
since, truncation points will affect the frequency response, 
they may be chosen on this basis. For a truncated sinusoid 
irrpvilse response, the frequency re^onse is given by. 


H(w) 



Cos (-^i-t), e3<p (-jwt) dt 


( 2 . 8 ) 


where, 

B = Tp/2 


( -Q. -w) 


si nXltw/O. ) Tip 
C^lhw) 


{ 2 , 9 ) 


where 


and 


-Cl = inpulse response sdLnusoid frequency 
p = number of cycles in iirpulse response 
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To eliminate the effects of d.c, offset in fault current it 
is required that, 

H(0) = 0 (2.10) 


which is 


2^ sin (jip) 
jfi " 


= o 


( 2 . 11 ) 


That is, the inpul ^ respon se is restricted to having an 
integral number of electrical cycles. The high frequency 
components are eliminated by analog low pass filter. 


2, 2. 2, 2 Selected Harmonic Filters 

While the differential equations of the line are being 
solved by nxmerical integration, the integration is carried 
out over a certain number of overlapping siibintearvals with 
required end points. This results in the elimination of 
certain harmonics and their multiples. 

An equation of the form, 

di 

v = R i t It -prr^ (2,12) 

X at 

is obtained for all fault conditions where the values of v. 
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found by appropriate choice and combination 
of sanpled voltage and current. In general, it is possible 
to calculate R and L, so that, any number of harmonics can be 
eliminated. For example to remove the third harmonic, above 
equation can be integrated once, over the interval O, c 
(a is an arbitrary constant) and over the interval Tt/S^ 
(a+Ti/l), the resulting equations, when added give : 


L 




a+iT/3 

4t/3 



i dt '1- 

X 


. cc-i.-ji./3 


'71/3 


i dt? = 

X 



(2.13) 


This equation is solved by numerical metliods in which the 
third harmonic and its mxiltiples are completely filtered. 


In general, to remove tv/o harmonics of order m and n 
and multiples thereof, the resulting equation will be ; 


2Ti/m 


27r/m+7t/n 


di, + 




1 


di. 


rc/Vf 


4" R 


27t/m 


I 


i dt + 

X 


2n/m-fn/n 

ir/n 


i dt 

X 


2r[/m 


V dt 


27T/m+7i/n 


n/n 


V dtj 

J 


(2.14) 
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2.2,2, 3 Least Square Error Filters [13^22] 

A formal procedure for smoothing the data is to fit 
a polynomial of degree m to {2n 1) samples of data. If 

m < (2n-i-l), a least square fit can be obtained. The necessary- 
time derivatives are then estimated by differenciating the 
polynomial , 

A cubic fit over seven points or a quadratic fit 
over five points or a straight line fit over three poin-ts are 
perhaps necessary to satisfy the redundancy requirement. If 
a waveform contains a decaying DC offset, the fundamental and 
a desired number of harmonic conponents are assimed and the 
least square error criterion is applied to determine the 
unknown parameters of the fundamental conponents. 


2, 2, 2, 4 Orthogonal Transform Filters 


Of these, the Fourier Transform filter is the most 
widely used, it suppresses all the unv/anted components and 
therefore, offers the best accuracy, but, it requires a data 
window of one full fundamental period [2,10 ], The orthogonal 
functions may include -the following : 
i) sine and cosine function - 


The Foufier transform F^(s) of cos n x is 


F^(S) - 


cos Tt X e 


4.T 


-j27TS X 


dx 


(2.15) 
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Since, cos 71 x is an even function 


F (s) = 


£ 


cosn X cos 2 ns dx 


(2.16) 


Using trignometric identies and solving, we get. 


Fc(s) « ^ sine sine 2 (s-^) 


(2.17) 


where, we define. 


sines = 


s xnK s 
ns 


(2.18) 


The Fourier transform F (s) of sinnx is given by 

3 




P^(s) = 


• ^ - j 2 tc sx , 

sinrtx e dx 




(2.19) 


which can be e2<pressed as. 


F^(s) = ^ sine 2 (s-|-) - ^ sin c 2 (s-t~) 


( 2 . 20 ) 


Sxaminaticn shows, .that, cosine/sine functions have a high 
rejection to dc and to high order harmonics. 


ii) Even and odd square wave - 

The iirportance of this function in computer relaying 
is, that, real and Imaginary corrponents of the physical 
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qiuantity sanpled are obtained trough addition rather than 
multiplication. 

The Fourier transform, even square wave 

is. 


Fj,(s) e 2 sines - sinc2s (2,21) 

The Fourier transform, ^’qCs) o£ odd square wave is, 

F^(s) *= -2j sincs.sinTTs (2,22) 

Examination shows, that, the square wave functions have a 
high rejection to d.c. and tend to limit high order harmonics. 

The properties of Walsh functions and Walsh trans- 
forraation [31 ] are used to extract the peak values of funda- 
mental frequency components of the voltages and the currents 
and then the inpedance seen by the relay is coirouted. 

Filters using sine/cosine functions, eitploying a data 
window of half cycle and less than half a cycle with tolerable 
error, have also been proposed [13,15], 

2, 2, 2. 5 Finite Transform Filters [23] 

The operational relationship given in eqn. (2.23), 
applies irrespective of the waveforms of voltage and current 
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for a solid faxilt in a transmission line : 


Vj, = z (P) ij, 


(2.23) 


where/ is the instantaneous value of relaying voltage 

i^ is the instantaneous value of relaying current 
Z(p) is the operational value of fault loop inpedance 

The voltage can be considered as the summation of an arbi- 
trary number of components v ^2 •••/ which are chosra/ 

so that/ at all times = "'^rl'^'^rZ"^ '** * Superposition can 
be used to relate the total current i^ to the components of 
current which would flow in the circuit in response to the 
separate application of each of the voltage conponents; i«e, 
if the currents that flow in response to the separate applica- 
tion of the voltage conponents v^^/ ^^^ 2 * '** ^rl' ^r2' '** * 

respectively/ the total current is given by 


^rl ^r2 


•i" m » m m 


(2.24) 


Each pair of conponents is related by an equation of the form 
of eqn.(2,23) / and/ considering the second corrponent, the 

time domain relationship of eqn. (2.25) applies : 

r2 


v^2 ~ ^ 


(2,25) 



20 


It is possible to transform eqn, (2,25) from time domain to 
the frequency domain by means of the Fourier transform inte- 
gral. 

The transformation, of eqn.(2,25) into frequency domain^ 
therefore/ results in the relationship given in eqn. 


Vj.2(jw)' = Z{jw) i^2 Uw) 


(2.26) 


where# 


''r2 (jw) = I v ^2 (-jwt) dt 

-(30 


1^2 (Jw) 


■r 

'i-CO 


i^2 (-Jwt) dt 


(2,27) 


(2,28) 


If# conponents at a particular extracted frequency are evaluated 
by means of the Fourier integrals, the fault loop iitpedance at 
the particular frequency chosen# is given by eqn, (2.29) 


If, for exanple, the component v ^2 equal to the relaying 
voltage v^ for all time inside a window T^ S. t <. T 2 


and 
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zero outsi^G^ thus window length being (T 2 -T^)/ 
given by, 

V^jOWg) = 7 (Jw^) 

2 

exp (-jwt) dt (S.S*^) 



The above eqmtion is referred to as the finite transform of 
the relaying voltage 

2. 2. 2. 6 Kalman Filters [24, 25] 

This metnod of filtering is especially suited to on- 
line digital processing because the noisy input (meas^arement) 
data is produced recursively. The non fundamental conponents 
in the current and voltage waveforms are designated as noise 
signals. As the noise signals depend upon the location and 
the type of faults, which are random in nature, the noise signal 
is considered a random process. The filter is initialized 
with an Initial estimate of signal and its error covariance. 
Then, as each measurement sample becomes available in real 
time, it is used to update or refine the filter's previous 


estimate 
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The noise signal in the voltage waveform is considered 
as a white noise seq[uence v/ith decreasing variance, A two 
state Kalman filter is used to extract the f undam 31 tal voltage 
phasOr, 


The noise signal in the current wavefojrm is considered 
as an exponential conponent plus a white noise sequence with 
decreasing variance, A three state Kalman filter is used to 
extract the fundamental current phasor. The third state 
accounts for the exponentially decreasing coirponent. These 
filters are expressed in the form of state equations. 

The error in this filtering process [24^25] has been 
reported to be less than 1% after half a cycle. The method 
appears to be suitable for high speed digital computer relaying 
based on post fault phasor (ie fundamental sinusoidal) quan- 
tities. 


We now give the overview of the digital distance 
relaying algorithm. 

2.3 • DIGITAL DISEANCS HSLAYIHG ALG0RITHJ4S 

These algorithms aim at extracting the fundamental 
power frequency corrponents of voltage and current from the 
complex post-fault waveforms and then, deteirmining the impedance 
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from these fimdemental Gon 5 >onents, A nirnber of algorithms 
have been proposed by several authors ['I “233 to 

calculate fault inpedance from the relay location to the point 
of fault. These algorithms can be grouped into four cate- 
gories which are described briefly below. 

2 ^ ^ Method 1 

In this method, the magnitude of inpedance is calcu- 
lated as the ratio of peak voltage to peak current and its 
argumait as the difference between the phase angles of the 
voltage and current, A few algorithms based on this method, 
as proposed by different authors [l~6] are given below, 

2. 3.1,1 Fourier Analysis Method - 

This algorithm was proposed by Ramamoorty [ 2 ], The 
sarrples for voltage and current are collected over one full 
cycle. The calculation of fundamental quantity is done as 
follows : 

f^ (t) ^ (2,31) 

^1 “ fQf2f^cosx-l-2f2Cos2x-r, . ,+f^cosnx) (2,32) 

b^ * ^(2fj^sinx+2f2Sin2x+...+2f^_^sin(n-l)x> 


(2.33) 
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where X is the sanpling interval (rad) givei^ by — 

XI 

XI is the nuiriber of intervals in one cycle/ 

^2 ••• sampled values of input signal 

over a period of one cycle. 

The factors cosx/ sinx ... cosnx/ sin(n-l)x are cons- 
tants and can be calculated and stored in the memory 
apriori as weighting functions on the sample values. 

The connroutation is done for both the voltage and the 
Gxarrent/ the voltage peak and current peak are calculated and 
then modulus of inpedance and its phase angle are evaluated, 

2. 3, 1,2 Predictive Calculation of Peak Values 

The peaks of voltage and cxrrrent are predicted before 
their actual occurance from the given samples [4] , 

Let/ V sinwt (2.34) 

Taking derivative of the equation (2.34)/ we get, 

V* ss w V T cos wt (2.35) 

pK 

From these two equations/ the peak value of sinusoidal voltage 
is given by. 



w 


(2.36) 
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If current is given by^ 


i = sin (wt-i-9) 


(2.37) 


Taking derivative of the eqn, (2,37?)* we get. 


i* = w Ipj^ cos (wt+<?) 


(2.38) 


The peak value of sinusoidal current is. 


^k 


/ . 2 / r \ 2 

IX + ( — ) 

w 


(2.39) 


Irtpedance modulus is given by. 


V 


peak 

^peak 


(2.4D) 


Dividing, eqn (2.34) by (2.35) 


V . sinwt 
vw _ pk 

v' ** Vpj^.coswt 


tan wt 


w 

V* 


^ -1 /VW% 

wt * tan 


(2.41) 
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Similarly for curren't sanple 

(wt-;-<?) = tan“^ (^) (2.42) 

Thus phase angle between voltage wavefoam and current wavefoim 
is^ 

= arc tan (^) - arc tan (^) (2.43) 

By using central difference theorem or baclcward difference 
theorem/ the irrpedance and its angle are evaluated. 

The predictation of the voltage peak and the current 
pealc using sanples and their first and second derivatives was 
proposed [?] / and also by using samples and sanpling interval 

[ S] • 

2.3,2 Method 2 

In this method/ the fundamental conponents of voltage 
and current are determined in the phasor form/ from which the 
real and imaginary parts R and X of the psuedo-iirpedance seen 
by the relay can easily be evaluated. A few algorithms/ based 
upon this method/ are described below. 
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2 » 3,2,1 Inpedance Calculation based on the Theory of 
Ozrthogonal Transform [lo] 

Since, inpedance is a cortplex function, the digital 
impedance relay determines the conplex value of the physical 
quantity, current or voltage, by sanpling with two orthogonal 
functions of the variable time. Two sets of such orthogonal 
functions viz, cosine/sine functions and even odd square wave 
fxanctions, have been discussed in section (2. 2, 2. 4), 

The real axis frequency spectrum is obtained by using 
the cosine function as a sanple function »i.e. real axis trans- 
form of physical quantity is multiplied to the Fourier trans- 
form F (e) of cosnx, similarly, the imaginary axis spectrum 
is obtained by using sine function as a sample function. 
Examination shows, that, these functions have a high rejection 
to dc and to high order harmonics. 

If the extraction of fiandamental frequency conponents 
is done by means of correlating the signal with the sine and 
cosine fxmction of the Fundamental frequency, and the data 
window is shorter than one cycle, the presence of a periodic 
Goxnponents in the signal give rise to large source of error. 

To minimize the error in calculating the fundamental corrponent 
for a data window equal to half cycle, the signals ought to 
be correlated with sine/cosine fxmction which have period equal 
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to data window length [16]. 

Let, Current and voltage input signals are, 

-(t/T ) -t/x 

l(t) — cos (w- 1— c) -j-I e -i- I e ^ cos(w t— r) (2,44) 

-I- a p p 

- (t/T ) -t/T 

V(t) = V, cos(w, t-a-{- )W^e ^ +V e ^ cos(w t-g) (2,45) 

X X a. p p 

First terms on the right hand side of equations (2,44) and 
(2,45) are steady state fundamental coupon ents. The second 
terms are well known aperiodic (that is dc) components decaying 
with time constant T_, The third term represents decaying 
oscillations induced by the fault. The real and imaginary 
parts of voltage and current are calculated for the angular 
frequency of W 2 corresponding to data window at which the 
spectrum of aperiodic conponent reaches minimum. It is known 
that. 


'1 


'd 

^d 



(2,46) 


(2.47) 


(2.48) 


where t]_ is the beginning of the data window Tw. 
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According to Fouxier Transform theory and correlating 

the signal with sine/cosine functions which have period equal 

to data windea?, Tw, that is we have, 

2 Tw ^ 


Y ~ .. iS . , 
d Tw 


P 

V ^ 

q Tw 


Tw 

r2 


Tw 
2 

Tw 

F 2 


Tw 

2 


V (t) cos W2'CdT 


V(t) sin W2'CdT 


(2.49) 


(2.50) 


I ^ ^ 

■^d Tw 


” Tw 


Tw 

r 2 


Tw 
" 2 

Tw 

r2 


■J Tw 
~ 2 


I(t) cos W2 TdT 


I(i:) sin W2 t dT 


(2.51) 


(2.52) 


where coefficients K and P are. 


K 


n (l - r ) 
r sin n r 


P = 


,^71(1 - r ) 
^injrr 


w^ 


and r = 


w,- 


Here, is the normal power frequency of 5 rsi^, and the 
angular frequency corresponding to data window. 
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Taking ' data window 

7T 

W^ 

written as/ 


half cycle, i.e. Tw 


equal to half the period, that is, 
, the coefficients IC and P can be 


K = 



( 2 , 53 ) 


Substituting the Values of coefficients K and P in the above 
equations, and using trapezium method of integration, the 
values of V,, V , I, and I are evaluated, 

d q-^ d q 

The impedance seen by the relay is obtained by dividing 
the psik of voltage by the peak of current, that is. 



-1- 

j V 

J , 




^d 

J 



v,i,-:-v I 
d . d q _q 

+ j 

V I. 

- V. 

T 

I,tl 
d q 


q 


Since Z = R+JX 
Therefore, 


V^I,-{-V I 
d d q q 

j t 


d q 


X « 


VI,- V,I 

. s. 

I, t I 
d q 


( 2 . 54 ) 


( 2 . 55 ) 
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2, 3, 2. 2 Irtpedance Calculation Using Least Square Fit [ 12 ] 

An example of a third degree polynomial fit to 
seven point data is given. The smoothed values of y^/ 
and Yq’* / for a given data set where, K = -3, -2, -1, 

0, 1, 2, 3, are given by. 




•2 ^ Y_3-^y3 ) +3 (y_2-^y2 ^ ^ 

2T 


(2.56) 




(2.57) 
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(2.58) 


These estimates are then used to determine the phasor represen- 
tation for {A 

2 2 2 ,r. cn\ 

y = yg -i- yg ( 2 . 59 ) 

where, y is real conponent, and y is imaginary corrponent of 
s c 

fundamental coirponent y. 


2 ,~,2 1 ,-,.2 
y cri (y) t — 2" (y / 
w 



w w 
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(yp -1- (yp (y^) cos (kwt) 

sjjr^ Oart) 


(2.60) 


and^ 


w(y ), 

(y‘) 


tr t < 

^ .-..T. ,X 

(wy* ) 




(y|) cos kwt - y^ 
(y|) sin kwt 


(2.61) 


where/ the subscripts 1 and 2 refer to estimates at sarrple 
times t^ and t 2 such that (tj^-t 2 ) = kwt/ k being and integer. 

An another algorithm/ based on the least error squares 
curve fitting technique/ assumes that input is composed of a 
fundamental frequency coitponent/ a decaying dc and harmonics 
coirponent of specified order [20 ] , The fundamental frequency 
component is determined with the above technique. The real 
and imaginary conponents of the voltage and the current are 
defined as 

^2v 

*2i 


« V 

cos 


sin 


p 


3v p 


V 


cos 

> *31 “ 

sin 

«1 
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The real and imaginary coxiponents of the impedance are given 
aS/ 


R 


Re (~) 


(X 


2v ^2i ^3v ^3i^ 


^^2i ^2i 


ic jc-'rr 

3i 3i 


X = Im 




^^3 v ^2 i ~ fgv 
^2i'^2i '*■ ^3i *^3p 


(2.61) 


(2.63) 


2.3.3 Method 3 

In tliis method, the line is modelled by differential 
equations, the numerical solution of which yields the value 
of R and X, The differential equation describing the trans- 
mission line (neglecting shunt parameters) is, 

V = Ri 4- L ^ (2.54) 

In the above representation, there is no need to suppress the 
D.C. offset as it is included in the representation of trans- 
mission line. A number of algorithm were proposed to solve 
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the ©ncjn. ( 2 , 64) numerically. In 1971^ Me Innes et al. [17] 
proposed an algorithm in which integration wa.s carried out 
over two successive time intervals. Integral equations are 
solved numerically using trapezoidal rule. The final e 3 <pre-, 
ssion for R and L are as shown below : 


R 




(2.65) 





( 2 . 66 ) 


where *v‘ and ‘i* are instantaneous value of the voltage and 
the current, *K‘ is the instant and ’h‘, the time interval. 

The above euqation assumes perfectly tra^nsposed transmission 
line and, it does not account for the shunt capacitance or 
series compensation. The fault resistance and the effect of 
power flow on the transmission line at the moment of a fault is 
also not considered. The data window is one cycle. 

In 1975 , Ranjbar and Cory [ 21] developed another 
technique, in which particular harmonic and their multiples 
are eliminated by taking appropriate integration intervals. 
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This technique has been presented in section (2,2,2). 

In 1979/ Smolinski [20 ] proposed an algorithm which 
accomodates both DC transient offset , transient high frequency 
components of line voltage and current signals without any 
additional filtering for these transients. The transmission 
line model used in this algorithm is a single equivalent tc 
section for the portion of the line from the relay location 
to the point of fault. The approach used in this algorithm 
is therefore one of including the sources of transient conpo- 
nents of the line voltage and the current signals in the 
irnpedance calculation algorithm/ rather than one of filtering 
the transient components out and, the, calculating the impedance 
of line from the remaining filtered signals. 


During post fault conditions, this single equivalent 
71 section model is described by equations (2,67) or (2.68) 


V 




a(i-ic) 
' dt 


(2,67) 


ic = ^ 
or 


dv 

dt 


V = i 


L _ R C ^ 

h at L dt 


- C 

JU 


d^v 


( 2 , 68 ) 
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By selBCting four successive sets of voltage and current 
sairples and replacing the derivatives in equation (2,68) 
by corresponding finite dif ferences/ one gets equation 
(2,69) from (2,68) for the four saiiple set. 



^K+1 " ^-1 

\ 

1 

tK-,-i-2v"K_i)n 





2h 

2h - 



■ 

'"k 


i ^-i 

K+2 K 

2h 

(yK+2-'"K) 

" 2h 

C^K+2~^ic-ri’^^l^ 
-T 1 





^K+2 

^K+3'"^K-1-1 

C''K-h3 - 

“ "2h"''' " 

C^K+3-^i<H-2"K+i 


1 



2h * 

h^ ! 


1 

] 

i 

^K-1-2 

^K-r3 


(yK+4-'"Kt2} 

" ■ 2h'''" 

h^ _j 



1 

1 

^K-t2 


(2,69) 


Ecfuation 2.69 can also be expressed in the partitioned form as 
shown in the ecruation (2,70), 


A 

B 


~p 


1 

i 

C 

D 


CP 


to 







L. ^ 


where/ 


- 

Lf 


J 


V 2 


II 


.^2 = 


— vr 


^Ktl 


>+ 3 _^ 


and 


(2.70) 
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A.,B .C 41) arc 2x2 stibniatrices 


sxncs, vector P and CP in the equation (2.70) are lin^rly 
related^ equation (2,70) can be reduced to equation ( 2 . 71 ) or 
(2.72) by matrix reduction teclinique^ and,. the unknovm line 
capacitance C is eliminated in the pr-ocess as well. 



(2.71) 


or 

- 1 1 
A p = V 


( 2 . 72 ) 


where. 


II 

1 A „ B C ' 

1 1 
]} V » 

y.1 ' 






— 


and P = 


R 



Erqoation (2.72) can be solved to give R^ and 
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2.3.4 M ethod 4 

In this method, the iirpedance is determined directly 
in phasor form from the ratio of the frequency domain values 
of the fundamental conponents of the voltage and current by 
finite transform method [23] , This method has been elabora- 
ted in the section (2. 2, 2, 5). The offline performance 
studies indicate, that acceptable accuracy can be coupled 
with an operating time of typically three quarters of a cycle 
of power frequency. 





CHAPTER 3 


PRQPQ:3BD dista nce RELAYING ALGO R ITHM 

3.1 INTRODUCTION 

In. this chapter the algorithm for the protection of 
the three phase transmission line using the digital corrputer 
is presented. The analog voltage and current signals are 
obtained through conventional P.Ts. and C,Ts. respectively^ 
and after passing through the appropriate circuitary^ these 
signals are t3mxiSHndh:^:.to the digital processor. The detection 
of the fault is done using cycle by cycle comparison of the 
voltage signals. After the detection of the fa\iLt, the appro- 
priate single phase relaying quantities are selected for sxibs- 
equent line irrpedance calculation. The voltage peak and the 
current peak of the sinusoidal waveforms and their relative 
phase difference are evaluated from the selected voltage and 
cuirrent signals and their first derivatives. The modulus and 
the argument of the inpedance are calculated from the division 
of the voltage peak by the current peak and the relative phase 
difference of the sinusoidal wavefoirms. Since/ the simulated 
fault data contain only fundamental frequency conponents/ no 
filters are errployed in this algorithm. The three zones 
protective schemes/ which are normally used for the protection 
of the transmission lineS/ are discussed. Different types of 
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irelay charac'tejri sties and the appropriate signals to obtain 
these characteristics in the pl="'e are discussed and these 
different threshold characteristics are shown in Tigs, (3.2a— 
3,3h) based upon phase comparison method. To avoid the 
possibility of mal— operation of the relay due to pov/er swing 
effects/ the relay characteristics such as quadrilateral/ 
elliptical and mho with blinders, which inclose tiie fault 
area coitpactly, are discussed. Finally/ the results, shovno in 
the tabulated forms, are examined and discussed, and the merits 
and demerits of the irtpedance algorithm are also examined. 

The digital relaying system consists of a general 
purpose high speed minicomputer or microprocessor coupled to 
a high speed A/D converter. In operation, sarrple and hold 
circuits internal to the multiplejcer, sairple phase to neutral 
voltages and phase currents 4Q times per cycle, i.e. in the 
interval of 0,5 ms. The A/D converter converts each analog 
signal to a digital equivalent whereupon its output is trans- 
mitted into memory via a direct access channel. 

The proposed algorithm is tested on DSC 1090 at 
I.I.T. Kanpur for all possible faults at several locations in 
the transmission line and also all possible protective relaying 
schemes are realized. 
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3.2 THEORY 

The proposed algorithm for the protection of the 
transmission line can be sub-divided into the follov/ing 
sequence of operations : 

1. detection of the faultj. 

2. classification of the fault, 

3. selection of single phase relaying quantities, 

4. dalculation of the irrpedance, and 

5. realization of any desired relay characteristic. 

Fault detection is performed using the cycle by cycle 
coirparison scheme as proposed by Mann and Morrison [ 5 ] • 

The counters such as NOFR, NOFY and NOFB are assigned to red, 
yellow and blue phase voltages respectively. The voltage 
sanples in each of the phases are compared to the corresponding 
samples in the previous cycle by subroutine COMP. If the 
difference between these two values for a particular phase 
exceeds a tolerance TOLV (i.e, 0,05 p.u.), the counter for 
that particular phase is incremented by one. When any of the 
counters exceeds a value CV (i.e, 5), the system enters into 
the fault detection mode. On the othorhand, if the difference 
between two voltage samples, found by above comparison is less 
than TOLV, the covin ter for that particular phase is decremented 
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by one, if it is not already zero. 

Zero sequence values of the voltage and the current 
are scanned by sixbroutine SAMP and are stored in the meinory. 

Once the fault is detected, the next job is to 
distinguish between two groups of faults, i.e, the phase 
faults and the ground faxilts. This is performed by subroutine 
FTYPS, in which cycle by cycle cortparison is done for zero 
sequence voltages. If the difference between this value and 
the corresponding value in the previous cycle exceeds the 
preset value TOLZ (i.e. 0,605p.u,), the fourth 
NOFG, assigned for ground, is incremented by one. On the 
O'tiierhand if this difference is less than the above preset 
value, GOvinter is decremented by one, if it is not already 
zero. When the counter value equals or exceeds to CZ (i.e. 5), 
the fault is detected as ground fault. 

After the classification of fault, the appropriate 
quantities of the voltages and the currents, i.e, the equivalent 
single phase relaying quantities are selected. The delta 
quantities for phase faults and phase voltages and zero sequence 
compensated currents for ground faults are needed for 

inpedance calculation. This is performed by subroutine SELCT, 
The aim of this subroutine is to determine which two of the red. 
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yellow or blue phases or ground are to be used to derive the 
equivalent single phase quantities. If a counter has value 5j. 
then the phase it represents is considered to be involved 
in the fault, if a counter has value 0 or 1, the phase is not 
involved in the fault, if the value is 2 or 3, it is undeter- 
mined. As phase faults include three phase symmetrical faults, 
line to line faults, line-line to ground faults, this algorithm 
treats them as the same. For exanple a RYG fault can be treated 
as RY or YG or RG fault, but it is preferable to treat it as RY 
fault rather than RG or YG fault since impedance calculation is 
more accurate for phase faults. 

The voltage peak and the current peak of selected sinu- 
soidal waveforms and also their relative phase difference are 
evaluated from the sarrples and their first derivatives, solved 
numerically using backward difference method (shown in Appendix 
All) , The modulus and the argument of the impedance are cal- 
culated from thr&~d4yisioii of voltage peak by the current 
peak and the relative phase di;^ference of the sinusoidal wave- 
forms of voltage and current respectively. 

The tripping decision is taken fro a the relay criteria, 
i.e. the relay gives a trip signal if the impedance lies within 
the tripping area of that relay for specified zone. 
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The time delays provided, for zone2 and zone 3 operations, 
are 10 ms and 20 ms respectively. Any desired relay charac- 
teristic is realized depending upon the requirements. 

The flow charts for the proposed relaying schemes 
are shown in Fig. (3.1a) - (3.1e). 

3,3 THREE ZONE PROTECTION OF THE TRAl'BMISSION LINS 

Different types of possible faults in a 3-phass trans- 
mission line are as follows ; 

1, three phase symmetrical fault, i.e, only one 
in no., 

2, single line to ground faxilts, i.e, 3 in no., 

3. line to line faults, i.e. 3 in no. and 

4. double line to ground faults, i.e. 3 in no.. 

Therefore, there are total of ten faults. These 
faults can further be classified into two categories such as ; 

1, ground faults - which include only single line 
to ground faults, and 

2. phase faults - which include rest of the faults. 

In general, the fault takes place involving more than one phase 
through fault resistance, which include arcing resistance of 
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the arc path^ tower footing resistance and resistance of ground 
return as the case may be, thus, the relay sees the impedance 
more than the actual value and hence over-reaching of the relay 
may take place. To avoid any maloperation of the relay because 
of the over reaching, normally the reach of the relay is set 
to 80%-90% length of the transmission line is being protected. 
Usually, a 3 zone protective scheme is adopted for the protec- 
tion of the transmission line, the first zone covers only upto 
80'/a~90yo length of the first section and operates instantaneously. 
The zone 2 covers entire laigth of first section and vj^to 30:4 - 
50% length of the second section and operates with some time 
delay and the zone 3 covers the entire length of the first and 
the second sections and operates v/ith some more time delay. 

The zone 2 provides the primary protection to the first section 
from 00% - 90% to 100% of its length, viz the laagth of the 
first section not being covered by the first zone, and for 
other parts of sections it provides back up protection. The 
zone 3 provides back protection to the entire lengths of 
the first and second sections and also acts as a starting ele- 
ment of the relay* This scheme is shown in Fig, 3,2, 

The voltage and the current sarrples at the relay 
locations are obtained through PTs and CTs respectively. The 
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replica inpedance or set inpedance is given by, 

-7 Cc T. r atio 

** 1 ^ P.T, ratio 


where is the line inpedance. The set irrpedances for all 
the three zones Z ^2 ^R3 respectively, are. 


Z 


R1 


0,8 X X 


C,T,__r^io 
P.T. ratio ' 


= (Z,^ t 0,5 X Z^2) ^ 


C,T, rat io 
F,T, "ratio 


^R3 




C^._j:at^ 

P.tT ratio 


where, and Z^2 inpedances of the first section and 

the second section respectively, 

Por Ijipalance calculation, the delta quantities of the 
voltages and the currents are used for the phase faults and 
phaas iTOltogos aisdczma .sequraoeUKSiireirsatei q^aseoouxieats are 

uaaad dtor'/taaa ground faults, viz- 



for phase fault 
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V 


Z = -. 


-1' K.I^ 
a o 


for grc’." .d faxill 


v/here is zero sequence current and K. is knoi-m as coiipen- 
sating current coefficient,- waich is given by^ 


K 


(Zo - Zl) 


z. 


where is the zero sequence irtpedance of the line and 

t 

is the positive sequence irtpedance of the line. 


3,4 TYPES OP DIST-nLTCE RELAYS 

Different types of classical relays and signals to 
obtain these characteristics are given in the following table, 

- "'Siplitude C^rtparison Phase Corrparison 


Characteristics 

Operate 

(S^) 


IZj^W 

Directional 

(Z-^+Z) 


Plain Irtpedance 


Restrain Measuring Polari- 
(S-) (Si) lizing 

„ __ _ ^ Csp ^ 


V_IZ^ 


V 

(Z-Zj,) 

(Zj,) 

. (Z) 

V 

IZ^-V 

iZr-;-v 


(z^^-z) 

(Zr+Z) 

.yi 'f* ’■> • 

87 ^ 

p 

8 / 


' ‘ 
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•bable co n'tln.ued 


Characteristics 


Angle Inpedance 
(ohm) 


Angle admittance 
(mho) 


Offset mho 


Operate 

(Si) 

Cotrparlson 

?.e strain 

(s^) 

Pte se Comparison 
Measu-fin^ pblarili- 
(S.;) zing (S^) 

2IZ^-V 

V 

IZ^-V 


(2Zj^-Z) 

(Z) 

(2 -Z) 

U. V 

'^R> 


2V-IZ^ 

IZ_-V 

R 

V 

(^R> 

<2Z-Zj^) 

(Zj,-Z) 

(Z) 



IZ^-V 

V-IZ 

o 

(Zr-^o) 

(2Z. ) 

(Z^-Z) 

(Z-Z^) 


As ainplitude conparator/ the tripping criteria is given by, 

I 

and at threshold. 



where S- and S„ are two signals, i.e, the operating and the 

X o 

re stirs, in ing to the relciy respectively 
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As phase conparator/ the tripping criteria is given by, 

o o 

~ 90 ^ * 

and at the threshold, f + 90^ ^ 

where 'I' is the angle between two signals S£ and to the 

rel^ ( f = - /Sp 

To obtain the reactance characteristic from the ohm 
relay characteristic, the angle of the polarizing signal is 
set to 90° i,e., the replica impedance angle is set equal to 
90°« Basically, this relay is used for ground faults as it 
is not affected by the fault resistance. This relay is always 
acconpanied by fault detecting and directional relays such as 
mho relay, which combines the function of starting and direc- 
tional feature in a single unit. 

To obtain the restricted characteristics and hence an 
another group of threshold characteristics, the threshold angle 
is set to an angle g, which is less than 90°, The conditions 
for the restricted threshold characteristics (phase comparison) 
are given by, 

-g < w < g for tripping. 
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and at threshold^ ^ = +3/ 

whera ^ <90'^ 


The another group of relay characteristics, obtained 
by restricting the classical ral^ characteristics, includes 
the following s 


1, restricted directioral, 

2 , restricted ohm, 

3, restricted mho, i,e. elliptical characteristic, 

4, Combination of restricted directional and restric- 
ted ohm, i,e, quadrilateral characteristic and 

5, mho relay with blinders. 


The threshold characteristics, as mentioned above, are 
shown in Figs, (3,3a) — (3,3'k) for the set irtpedance For 

the three zone protective scheae, Z is replaced by appiopriate 
set inpedances, i.e. it is replaced by and Z^^ for 

first, second and third zone respectively. 


Flow charts for 3 zones protective schemes, such as 
plain inpedance with directional unit, mho relay with offset 
characteristic as zone 3 unit, reactance relay alongiirith mho 
relay as directional unit, elliptical relay, q^iadrilateral relay 
and mho relay with blinders, are shown in Figs. (3.4a) -(3. 4f) , 




FIG. 3.2 THREE ZONES PROTECTIVE SCHMES 




SIGNALS 

51- ^IZr-V 

52- ^1Zf^*V 

AT THRESHOLD 90® 
OTHERWISE -90® < f < 90® 


FIG.3,3b PLAIN IMPEDENCE RELAY 




r IG. 3.2.C 


OHM RELAY {U. angle 
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/ y\ 

( 

j /xmpj 


y 


Tr^ 


impedance relay ) 

SIGNALS 

s\ “V 

S2— V 

AT THRESHOLD Y** 90® 
OTHERWISE -90®<‘f<90® 


IG. 3.2.d MHO RELAY (i.e. angle admittance relay ) 



SIGNALS 
Si'-IZr-V 
S'2-KIZr+V (Ktt ) 

AT THRESHOLD <p«t90® 
OTHERWISE -90®<<f<90® 


FIG. 3.241 OFFSET MHO RELAY 
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SIGNALS 
Si—IZr -V 
Si — IZf, Zso® 

AT THRESHOLD 

tp *1 90 

OTHERWISE - 90® < <p< 90® 



SIGNALS 

Si^IZr 

S2-i-V 

AT THRESHOLD 

and 0<90® 
OTHERWISE 


FIG. 3.2.g RESTRICTED DIRECTIONAL RELAY 




SIGNALS 

S'— IZr 
AT THRESHOLD 
tp*!fiand 0<90® 
OTHERWISE 


r* i 


-i-jK RpciTRICTED OHM RELAY 



SB 



FlG.3.2.i ELLIPTICAL RELAY (Restricted MHO relay) 


& IZj, SIGNALS 

wj S'l — IZR 

S2— izr-v 

^2 ^/m THRESHOLD fllneKSO* 

P"^— OTHRWISE ^ 

„,«) 0U«R«TE«AL BEL.«R..««a .hmE *•.««« 



SIGNALS 

sJMIZr-V) 

S'2-V 

S'3--IRr/ -o£ 

-V) 


FtG 3 2M 


mho with 8UNCCRS 
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FIG. 3.4.0 PLAIN IMPEDENCE RELAY 
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CALCULATE ANGLE, TH1, BETW- 
EEN SIGNALS,(IZf^l-V)AND V 
CHOOSE L B, AND CONSTANT k 



CALCULATE ANGLE ,TH2,BETW- 
EEN SIGNALS, 


(IZro-V )ANpV 



FIG 3.4.d ELLIPTICAL RELAY 
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FIG.3.44 MHO RELAY WITH BUNDERS 
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3.5 RSdULTS hND discussions 

The proposed three phase relaying scheme has been 
tested for the transmission line (App, I) with the simulated 
fault data for all types of faiilts.. i.e. single line to 
ground faxilts, double line to ground faults, line to line 
faults and three phase symmetrical faults. The errors in the 
calculation of modulus of the inpadance and its phase angle 
are approximately within 10% and t 12 deg, respecitvely. 

The calcuLated line impedances and actual impedances, for 
several locations of faults are shox-m in tableS-l, 

In the simulated fault data, the fault impedance is 
neglected and . dc offset and h. f , transients in the post 
fault waveforms are also neglected, thus no filter is being 
realized. 

After the occurance of the fault, only seven samples 
are needed for fault detection and subsequent impedance 
calculation. The tripping signal is obtained based xpon the 
criteria of any desired relay characteristics, several three 
zones relaying schemes have been realized. The time delays 
provided, for zone 2 and zone 3 tripping, are 10 ms and 20 jns 
respectively. The operating time of the relay is found to be 
approximately 10 ms for zone 1 tripping. 
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CHAPTER 4 


CPN CLU SIOI-jS 


This chapter is aimed at reviewing the siginif leant 
results obtained during the course of this work and making 
a few suggestions for future scope in this field of prote- 
ction. 


In chapter one/ it is stated that larger amovint of 
power can be transmitted without loss of stability/ when the 
fault is cleared in a shorter time. Thus one of our major 
aims is to make the relay fast in operation alongwith other 
merits. The digital coirputer protection of transmission 
lines possesses several merits such as faster in operation, 
flexibility in obtaining any threshold characteristics and 
self checldLng facilities against hardware failures. In the 
same chapter, the overview of the wozks done by several 
authors in the field of digital protection so far, have been 
outlined in brief. 

The filters are enployed to extract the fundamental 
frequency components from post fault conplex waveforms as the 
distance inpedance relay algorithms are based upon only the 
fundamental frequency components of voltage and exarrent, some 
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of these algorithms and filters are presented in chapter 2 , 

The algorithm for the protection of three phase 
transmission lino based upon the predictive calcxilation of 
the jrrpedance is being discussed in detail and also tasted 
for all tj^es of faults on transmission line at different 
locations. The line irtpedance calculation yields a high 
accuracy, i,e. within t 10%, This :*.s inspite of considering 
only the first terra in the expression for numerical differen- 
tiation, Even, better accuracies can be realized by consi- 
der 5.ng a few more terms in the above expression. A number 
of 3 zones protective relaying scher.'.es such as plain irpe- 
dance relay with directional unit, mho relay with offset as 
3rd zona unit, reactance relay with rriho as starting element, 
quadrilateral relay (i,e, restricted directional and restricted 
ohm relay), elliptical relay (i.e, restricted mho relay) and 
mho relay with blinders, have been realized with time delays 
being 10 ms and 20 ms for zone 2 and zone 3 tripping. The 
special relay characteristics such as quadrilateral and 
elliptical, which are suitable for long heavily loaded line, 
offer the best threshold characteristics as they enclose fault 
area corrpactly. 

The operating time of the rel^, for zone 1, is 
found to be approximately 10 ms for different types of faults. 
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Tile actual operating time of the relay may include the time 
delays caused in A/D converter^ multiplexer and appropriate 
filtering. 

This algorithm for inpedance calculation takes 
lesser time than that in other algorithms as the peaks of 
the sinusoidal voltage. and current are predicted before 
their actual occurance, but this algorithm accepts only 
fundamental conponents of the voltage and the current signals. 
Although, theories of filters have been discussed but, in 
the present work, no such, filters are enployed as the simu- 
lated. fa^alt data contain only fundamental components. 

The future work may include the harmonics cortponents 
in addition to the fvuidamental frequency corrponsnts in the 
simulated fault data and realization of appropriate analog 
and/or digital filters to extract the fundamental frequency 
components. The work may also be extended with the inple- 
mentation of actual scheme by designing conplete hardware so 
as obtain digital on line protection* 

The installation of single unit digital corrputer in 
power station finds the applications in the fields, besides 
the protection, of the transmission lines, such as protection 
of the generator against the loss of excitation and windings 
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faults/ control of reactive power (i,e, control of terminal 
voltage) and protection of transformers- bus-bars and other 
equipments installed in the sub— station and. also load dispatching 
and data logging. 

A very fast relay scheme can be developed using 
travelling wave phenomenon wherefrom the post fault waveforms 
of the voltage and current fundamental conponents are 
suppressed by appropriate filters thus working v/ith the 
transient components. 
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APPENDIX I 


Transmission line (Singraxjli to Kanpur/ 400 ■ kV SOj: 
data, 'Tarhich -"is 'Used fon testing the_pix3posed protoctive 
schema . 

Base : lOO KVA 

Length of transmission line : 456 Kms 

Line parameters : 

Positive sequence resistance/ = 0.00927 p.u. 

Positive sequence reactance/ = 0.09427 p.u. 

2aro sequence resistance/ = 0.08117 p.u, 

Zero sequence reactance/ =* 0.27419 p.u. 

Source inpedance (3jc2l0 MVA) : 

Zero sequence reactance/ ** 0,016 

Direct axis s\ibtransient reactance/ X^‘ = 0.03 p.u. 

Direct axis transient reactance/ = 0,0412 p.u. 
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APPENDIX II 


DIFIERSNIIaTIONS formulas 


The basic central difference suppression for the 
derivative 


H yjc 


(-All.l) 


where # jx and 6 are the standard notations for the opera - 

encing 

tions of baclcward differ-/ / central differencing and averag- 
ing respectively and h is step size, and y stands for v or i. 

Using the first term only 


*> yi =■ 2 <y)c+i - y^-i’ 


(All. 2) 


and with a second-term 


yi = - n ^yk+2 I y^+i 


2 1 
3 ^k-1 ^ y- 


12 -^k-2 


) (All. 3) 


For baclcward differences 


= (V** I v^- 


(All. 4) 


Using the first term only. 


*"yk = yjc - yk-i 


(All. 5) 
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and with a second term 

I ^^k " ^k-2^ (All. 6) 

Note that (All, 2) and (All. 3) for involve sanple values 
at times later than whereas (All. 5) and (All, 6) do not. 
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